clc; clear;

% 固定点的特征值稳定性分析

ksc_val = linspace(-1.4, 1.4, 500);   % F1 的 k*sin(c) 的范围
mu_1  = linspace(-1.4, 1.4, 500);    % mu 的范围
mu_2  = linspace(0.4, 3.6, 500);    % mu 的范围

% 初始化存储矩阵：1 表示稳定 (CST)，0 表示不稳定 (UST)
stab_F1 = zeros(length(mu_1), length(ksc_val)); 
stab_F2 = zeros(length(mu_2), length(ksc_val));


% 遍历参数网格并计算稳定性
for i = 1:length(ksc_val)
    for j = 1:length(mu_1)
        mu1 = mu_1(j);
        ksc = ksc_val(i);

        % 对 F1 的特征值
        lambda2 = mu1;
        lambda3 = mu1 + 2*ksc;
        if (abs(lambda2) < 1) && (abs(lambda3) < 1)
            stab_F1(j, i) = 1;
        end
    end
        
    % 对 F2 的特征值
    for j = 1:length(mu_2)
        mu2 = mu_2(j);
        ksc = ksc_val(i);
        lambda2 = 2 - mu2;
        lambda3 = (2 - mu2) + 2*ksc;
        if (abs(lambda2) < 1) && (abs(lambda3) < 1)
            stab_F2(j, i) = 1;
        end
    end
end

% 绘制图像
figure;
subplot(1,2,1);
imagesc(ksc_val, mu_1, stab_F1);
set(gca, 'YDir', 'normal');
colormap([1 0 1; 0 1 1]);
title('F1 的稳定性');
xlabel('k*sin(c)');
ylabel('\mu');

subplot(1,2,2);
imagesc(ksc_val, mu_2, stab_F2);
set(gca, 'YDir', 'normal');
colormap([1 0 1; 0 1 1]);
title('F2 的稳定性');
xlabel('k*sin(c)');
ylabel('\mu');
